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_Thermo-Magnetic Study of Ferromagnetic Minerals 
Contained in Igneous Rocks* 


By Syun-iti AKIMOTO 
(Geophysical Institute, Tokyo University) 


Abstract 


The change with temperature in intensity of magnetization 
of ferromagnetic mineral grains separated from the igneous rocks 
was measured in a strong magnetic field by means of a magnetic 
balance. The intensity of saturation magnetization and the Curie 
point of ferromagnetic minerals of titanomagnetite series, which are 
commonly found in natural igneous rocks as the solid solution of 
TiFe,O, and Fe,;0O:, were examined in taking their chemical and 
crystallographic properties into consideration. Linear relations be- 
tween each two of the chemical composition, the lattice parameter, 
the Curie point and the intensity of saturation magnetization were 
found as the first approximation. Besides, it was found that there 
exist natural ferromagnetic minerals in ilmenite—hematite series 
having the peculiar type of their thermo-magnetic curve. 

§ 1. Introduction. 

The magnetic properties of rocks and minerals have been studied by a number 
of investigators, especially in connection with the local geomagnetic anomalies and 
their variation with time.“ In recent years, the phenomenon of the magnetic 
polarization of rocks im situ has attracted the attention of many geophysicists, since 
it may .be used as a certain kind of fossil which can represent the ancient geomagne- 
tic field. In other words, the magnetization of rocks could be used in order to find 
the secular variation in the geomagnetic field throughout the long geologic time. 

It will be obvious that the magnetic properties of rocks are almost due to the 
magnetic properties of the ferromagnetic minerals contained in them. Hence it 
seems that the most fundamental way of clarifying the magnetic properties of rocks, 
on the basis of physics of matters, is to examine the magnetic properties of ferro- 
magnetic minerals by taking their crystallographic and chemical properties into 
consideration. If this could be done, various magnetic properties of rocks could be 
reduced to those of the rock-forming ferromagnetic minerals and then the former 
would be described fairly systematically. 

The studies of the magnetic properties of rocks made so far have mainly dealt 
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with those in a weak magnetic field such as the geomagnetic field, since the result 
of such studies may be directly connected to the magnetic behaviours of. rocks zu situ. 
However, such magnetic properties in a weak magnetic field are the Btructines sensi- 
tive properties, upon which small stress or grain size of the ferromagnetic minerals 
have rather large influence. 

At the present stage of knowledge of magnetism of rocks, it seems that the 
first thing to be done in order to establish a general catalogue of rock-magnetism 
will be to find the general magnetic properties of ferromagnetic minerals where 
chemical composition and crystal structure are given. For this purpose, the structure- 
insensitive magnetic properties such as the saturation magnetization and the Curie 
point, which are mostly independent of internal stress and grain size, efc., and 
dependent only on chemical composition and crystal structure, will represent the 
primary magnetic character of the rock-forming ferromagnetic minerals. 

Several studies along this line have already been reported.“ Although they 
are very useful, they are still insufficient for getting a general view of the funda- 
mental properties of rock-magnetism. In the present work, the author attempted to 
describle as systematically as possible the structure-insensitive magnetic properties 
of the ferromagnetic minerals separated from a number of typical igneous rocks in 
relation to their chemical composition and crystal structure. 

.It is a well-known fact that the ferromagnetic minerals in rocks are the metallic 
oxides which are chiefly composed of FeO, Fe,0; and TiO... The magneto-chemical 
study of these ternary system has already been performed in fair detail by E. Pouil- 
lard by examining artificially synthesized specimens.“ In the present study, the 
experimental results obtained from the natural ferromagnetic minerals are examined 
in reference to this result obtained by E. Pouillard. 

§ 2. Specimens. 

(A) Titanomagnetite series. 

The specimens examined here are 15 in total number and each of them is the 
ferromagnetic mineral separated from the typical igneous rocks, being called “ magne- 
tite” petrologically in the broad sense of the word. Their magnetic properties in a 
weak magnetic field, such as magnetic susceptibility, have already been described in 
the previous paper.“ In Table I, the chemical composition of these ferromagnetic 
minerals which were analysed by T. Katsura are summarized together with their 
locality, petrological description and crystallographic parameter.: Plotting the chemical 
composition of these specimens on the FeO-Fe,0;-TiO, ternary system shown in Fig. 
1, we find that they are distributed around the TiFe,O,-Fe,O, line. According to 
Pouillard, the equilibrium among FeO, Fe.O; and TiO, at the room temperature takes 
place at the points, D, E, F and G and the complete line EG and the limited ‘line 


between C and D in the diagram in Fig. 1. In this figure, the line EG represents 


the complete solid solution between TiFe,0, and Fe;O,; and the limited line CD 
represents the solid solution between TiFeO; and Fe,O;. The former. and the. latter 
will be called here the titanomagnetite series and the ilmenite-hematite ‘series 
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respectively. 

The X-ray analyses of these 
specimens were carried out by means 
of a wide-range X-ray spectrometer, 
“Norelco,” with Fe Ka radiation. The 
results are summarized also in Table 
I. It was found that the crystallo- 
graphic structure of all the specimens 
belongs to the spinel structure, the 
lattice parameter, a, varying from 
about 8.42A to 845A. 

From these chemical and crys- 


tallographic data, it may be concluded Fo Fad. FeO. 
Fig. 1. Chemical composition of ferromagnetic min- 
: A erals separated from igneous rocks, repre- 
of the minerals belonging to the sented on a FeO-Fe,0,-TiO. diagram in Mol 


titanomagnetite series; z.e. the solid percent. 
solution of magnetite with TiFe.O,. 
(B) ILlmenite-hematite series. 
It seems likely from the above result that the principal ferromagnetic minerals 


that these 15 specimens chiefly consist 


in rocks are generally the titanomagnetite. There are however several cases that a 
fair amount of minerals of ilmenite-hematite series are contained in rocks. The 
specimens examined here are six in total number; their chemical composition, petro- 
logical description and crystallographic characteristics are summarized in Table II, 
where No. 1’ specimens is the ilmenite ore minerals and No. 2’, No. 3’, No. 4’ specimens 
are the ilmenites, petrologically so-called, separated from the typical igneous rocks, 
while No. 5’ and No. 6’ specimens are the peculiar ferromagnetic minerals separated 


thermomagnetically from those rocks which show the phenomenon of the reverse 
thermo-remanent magnetism. 

Representing the results of the chemical analysis, which was performed on 
No. 2’, No. 3’, No. 4’ and No. 6’ specimens by T. Katsura and J. Osaka, as the points 
on the diagram of the ternary system shown in Fig. 1, it was found that these 


Table Il. Chemical composition and crystallographic and magnetic properties of the ferro- 


, Chemical composition* 
No. Locality Petrological desciption |" Fe,0,; FeO TiO: Fe.,0; FeO: TiO. 
| wt.% : mol.% 
i | Inariana, Gunma Ilmenite ore mineral 
2' Sinkamibasi, Kagosima Phenocryst ilmenite 
in dacite pamios 21.9 932.8 45.8 11.8 38.8 49.4 
3’ | Kaziya, Yugawara Phenocryst ilmenite 
ia dacits Gbnlaicn 16.1 39.9 44.0 8.3 46.0 45.7 
4’ | Matukizawa, Asio Phenocryst ilmenite . 
in dacite pitchstone 12.4 37.% 49.9 6.3. 42.8" 5029 
5’ | Matukizawa, Asio Phenocryst ferromag. 
min. in dacite pitchstone 
6’ | Hutatudake, Haruna Phenocryst ferroma 
’ g. 
min. in dacite pumice 43.15 24 ee Oene 26.9 33.2 39.9 
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specimens are distributed around the TiF eO;-Fe,0, line. The results of the X-ray 
analysis by “Norelco” with Fe target are given in Table II. It was found that these 
specimens have the rhombohedral structure and their crystallographic parameters, 
4, and @,,,, are distributed between those of ilmenite and hematite. 


§ 3. Apparatus for magnetic measurement. 

Since the content of ferromagnetic mineral in igneous rocks is usually less than 
several rercents, it is rather laborious to separate out a great amount of ferro- 
magnetic minerals from the mother igneous rocks. For the sake of practical conveni- 
€nce an apparatus which can be used with a sufficient accuracy for measuring the 
magnetic properties of a small number of specimens was especially deviced. The 
apparatus is the same in principle as the magnetic balance which is generally used 
for the measurement of the magnetic susceptibility of paramagnetic or diamagnetic 
substances. However, the mag- 
netic balance was improved in 
some points so as to be suit- 
able for the measurement of 
the ferromagnetic substance. 

The apparatus is illust- 
rated schematically in Fig. 2. 
The specimens, S, which is 
inserted in the fused silica 
bascket, is suspended by a 
quartz spring, Q, in the middle 
of the pole-pieces of the elec- 
tromagnet. When the inhomo- 


Fig. 2. Schematic view of measuring apparatus. 


geneous magnetic field was Qtr ae ok quartz spring 
produced by the electrontagnet, Siyscc wee te specimen 
_ see nitee « Ia adios BW cc Not 0c Pate electric furnace 
pt BOGE LOLEE.: Hs ~ a ION ee electro-magnet 
on the specimen is expressed Ty Tet: OOP ee thermo-junction 
by 
magnetic minerals of ilmenite-hematite series. 
Crystal parameters Curie point ip 
ia as “| at room 
apn ark H Cc temp. 
: ~ ; 7 
5.540A 54°44’ 100°C 150°C 0. 5 &-m.u. The numerical values were calculated 
ae in taking the presence of ‘some 
Cc / 
as aio = Nes amounts of hypersthene into consider- 
5.523 54°54’ 100 140 0.2 ation. 
' **H[ represents heating process. 
5.522 ae er a ade? C represents cooling process. 
5.483 55°02’ 230 mk 


5.480 55°08’ 250 | 20 
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2 0H 1 
F=m:o-—~~ (Tf) 
Oo = FT (2) 


where o, x. m and H are the specific magnetization of the specimen, the specific 
susceptibility of the specimen, the mass of the specimen and the magnetic field acting 
on the specimen respectively. The saturation magnetization of the specimen can be 
measured provided that the applied magnetic field is strong enough to saturate 
magnetization of the specimen. This magnetic force acting on the specimen is 
determined approximately by measuring the amount of the elongation of the length 
of the quartz spring with the aid of a microscope. It was ascertained that the error 
expected by the measurement of the elongation is less than 5% in the present study, 
provided that the elongation of the quartz spring is not so large, ze. about 1 mm. 
Hence it seems that this apparatus is sufficiently useful for the thermo-magnetic study 
of ferromagnetic minerals. 

Thus, the change of the magnetization of the svecimens with temperature in 
a strong magnetic field can be measured easily by putting the specimen into the 
electric furnace, shown in Fig. 2. In the present study, the measurement was made 
under the atmospheric pressure. The absolute intensity of magnetization of specimens 
could be determined by calibrating the instrumental constants.. 

The intensity of the magnetic field acting on the specimen in thermo-magnetic 
study is 2000~30000e. for all specimens. When the demagnetizing factor owing to 
the shape of the scecimen is taking into consideration, this intensity of the field will 
be sufficient for saturating the magnetization of specimens. Hence, the thermo- 
magnetic curve measured by the present apparatus may be considered to represent 
approximately the change of the spontaneous magnetization with temperature. 

Since the spncimen is very small, the uniformity of the temperature within it 
will be almost perfect, but, an error owing to discrepancy of temperature of the 
sample from that of the thermo-junction amounts to a few degrees. This error may 
be within the limit allowable for the present purpose of study. 

As mentioned above, this simple apparatus could be used for the study of the 
thermomagnetic analysis of ferromagnetic mineral with a fair accuracy; even for the 
magnetization of a single grain of the ferromagnetic mineral. On the other hand, 
an ordinary magnetic balance was used for measuring the fairly accurate intensity 


of magnetization of specimens at the atmospheric temperature by means of the null 
method. 


§ 4. Change in Magnetization with Temperature—Thermo-Magnetic Curves. 


The intensity of magnetization of all specimens described in § 2 was measured 
at continual temperatures from the atmospheric temperature up to the Curie 
point by means of the above-mentioned apparatus in the atmosphere during both the 
heating and the cooling processes. Since the thermo-magnetic curves of the specimens 
of the titanomagnetite series and those of the ilmenite-hematite series are appreciably 
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Cifferent from each other, it will be convenient to discuss about each of them 
separately. 

(A) Titanomagnetite series. 


The measurement was carried out for several grains which were sampled 
randomly from each specimen. The thermo-magnetic curve thus obtained in 
magnetic field of 2000~30000e. will be regarded approximately as the change of the 


i) 1 Hi 

a — Mt ma >—* Reheating process 

ok % =—— Recooling process 
™ 
05 Qs US 
0 0 0 — 4. re 
P ° 200 400 oo TC 0 200 400 wo TC 

, No. 1 H=17000e. No. 2 H=24000e. No. 3 H=26000e. 


of ne 


° 20 “~ a TC *s 20 400 mo TC 0 a) 0 oo TC 
No. 4 H=21000e. No. 5 H=33000e. No. 6 H=24000e. 
in 10 
) ‘ 05 
*o 20 400) mo TC °0 200 40 oo TC 
No. 7 H=30000e. No. 8 H=24000e. No. 9 H=17000e. 
44 
1p e, 
10 oF 
05 05 
0) 0 
Tn 400 wo TC 9 bd 400 wy TC 0 20 400 wo TC 


200 
No. 10 H=30000e. No. 11 H=21000e. No. 12 H=24000e. 


8 S. AKimoTo 


KK 
\ 


=} 


ce ae 
0 20) 400 oo TC 0 200 400 oo TC 


No. 13 H=28000e. No. 14 H=30000e. No. 15 H=17000e. 


Fig. 3. Thermc-magnetic curves cf ferromagnetic minerals of titanomagnetite series. 


spontaneous magnetization with temperature. 

The results of the experiments are shown in Fig. 3. As will be seen in these 
figures, the type of these thermo-magnetic curves“ is very similar to that of the 
ferrite of magnetite type. However, the thermo-magnetic curve of some specinmiens 
changes irreversibly with temperature in heating and cooling processes. This is the 

_same phenomenon that has been found previously in the measurement with temper- 
ature of the change of the specific susceptibility of these ferromagnetic minerals in 
a weak magnetic field in vacuum.“? The specimens of the phenocryst minerals in 


the original rocks were classified into “reversible type” and “irreversible type” with 
respect to the change of their spontaneous magnetization with temperature, just as 
in the case of their magnetic susceptibility in a weak magnetic field. The reversible 
and irreversible types in the former correspond to the reversible and irreversible 
types in the latter, except Ne. 7 specimen. It is noticeable for all specimens that the 
degree of the irreversibility is more remarkable in the magnetization in a weak 
magnetic field than that in such a strong field as the present study. The thermo- 
magnetic curve of the irreversible specimen in re-heating and re-cooling processes 
coinsides with that of the initial cooling curve, as shown in Fig. 3 for No. 3 and No. 
4 specimens. These facts may show that the initial: cooling curve is more stable 
under the present experimental condition. This has been noticed in the previous 
experiments on the thermal change of the magnetic susceptibility in a weak field in 
vacuum. In the above-mentioned irreversible type special the Curie point de- 
creases apparently by heat-treatment. Such specimens as No. 12, No. 13 and No. 15 
have somewhat different characteristics of the thermo-magnetic curves; the increase 
of the Curie point was found in the thermo-magnetic curve in the cooling process in 
comparison with that in the initial heating process. This seems to be due to some 
chemical reaction such as the oxidation owing to the smallness of the grains and the 
experimental procedure carried out in the atmosphere. 


It has been reported by Pouillard that the Curie point of titanomagnetite 
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series specimens decreases according as the content of TiFe.O, increases in the solid 
solution between TiFe,O, and Fe,0,. Hence, it may be inferred that the decrease of 
the Curie point of the irreversible type specimens is caused by the increase éf the 
content of TiFe,O, in the Fe;O,-TiFe,O, specimen owing to the heat treatment up to 
about 600°C. It will be found in Table I that the crystal parameter of the heat- 
treated specimen is a little larger than that of the fresh specimen of the same 
mineral. As will be discussed later, the amounts of the increase of the lattice para- 
meter is nearly in accord with the theoretical value derived from the amount of 
decrease of the Curie point, where the rate of increase in the lattice parameter 
corresponding to the decrease in the Curie point is estimated to be about 0.015 A/ 
100°C. 

(B) Lmenite-hematite series. 

The thermo-magnetic curves of the specimens belonging to this series are 
summarized in Fig. 4. As will be seen in these curves, some of the specimens 
changes stepwise in the heating process. This stepwise changes of such specimens 
as No. 1’, No. 2’, No. 3’ and No. 4’ will suggest that these specimens contain a certain 
amount of titanomagnetite. An important characteristic of all svecimens of this 
series will be the peculiar type of the thermo-magnetic curve, that is their magneti- 
zation changes with temperature almost linearly over a fairly wide temperature range. 
It was ascertained especially for No. 6’ specimen that the magnetization changes 
linearly from the Curie point down to the liquid air temperature. 

It may be noticed in Fig. 4 that all specimens of this series have a ferromag- 
netic character though the intensity of some of their magnetization is not so strong 
as in the case of the ordinary ferromagnetic substances. It will be a remarkable 
fact that there exists a ferromagnetic phase in this ilmenite-hematite series which 
has been considered hitherto, even in Pouillard’s recent paper, to be nonmagnetic. 

On account of the heat treatmeat in the atmosphere, the thermo-magnetic 
curves of No. 1’, No. 3’ and No. 4’ specimens showed the irreversible changes. This 
change may be due to development of titanomagnetite phase in these specimens, 
which may result in an increase in the intensity of magnetization. 

§ 5. Saturation Magnetization and Curie Point. 

The intensity of saturation magnetization and the Curie point of all specimens 
will be discussed here, since they are the intrinsic magnetic properties of the speci- 
men which chiefly depend upon the general chemical composition and the crystal 
structure, being independent of the grain size and the internal stress. 

(A) Titanomagnetite series. 

The magnetization curve of the specimens of the titanomagnetite series was 
measured at the room temperature by means of an ordinary magnetic balance with 
null method. The intensity of saturation magnetization was determined by the 
following empirical formula by extrapolating to H=.o, 

o=0,( 1-4), (3) 
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Fig. 4. Thermo-magnetic curves of ferromagnetic minerals of ilmenite-hematite series. 
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where A is the constant determined experimentally for every specimen. The intensity 
of saturation magnetization tabulated in Table I varies from about 30 é€.m.u./gr. to 
about 90e.m.u./gr., being always less than 93 e.m.u./gr. which is expected to be the 
intensity of saturation magnetization of pure magnetite at the room temperature.“ 

On the other hand, the Curie point in each specimen was determined from the 
thermo-magnetic curve described in § 4. The Curie point values, tabulated also in 
Table I, spread over the temperature range from 300°C to 500°C. The Curie point, 
@, thus determined is in good agreement with the Curie point of the same specimen, 
0’, which was determined from the susceptibility in a weak magnetic field. 

Since the titanomagnetite is generally considered to be a solid solution between 
TiFe.O, and Fe;O;, it may be expected that there exists a close relation among the 
chemical composition, the lattice parameter, the Curie point and the intensity of 
saturation magnetization. The relation between the chemical composition and the 
lattice parameter is shown in Fig. 5, where the content of TiO, in Mol percent was 
expressed as the mezsure of the substitution of Fe*? ion by Fet? ion and Ti*4 ion in 
titanomagnetite. The general tendency that the lattice parameter becomes larger as 
the content of TiO. increases is found in this figure. However, the relation between 
TiO, and the lattice parameter for the natural titanomagnetite grain is a little less 
definite in comparison with that for the artificially synthesized specimen as obtained 
by Pouillard. Since the chemical composition of these natural titanomagnetites are 
deviated from the exact Fe,O,-TiFe.O; series, as illustrated in Fig. 1, it seems likely 
that the discrepancy mentiond above is due to the condition that some amounts of 
TiO, and FeO in the natural minerals 
make the phase of ilmenite which 
intergrows into the phase of titano- 
magnetite. 

Since both TiFe,O, and Fe,0, 
have the inverse spinel structures, the 
structure of any titanomagnetite may 
be generally expressed by 

aFe**(Tit*Fe**)O,(1—2) Fig. 5. Lattice constant of ferromagnetic minerals of 
Fe**(Fet*Fe**)O, titanomagnetite series as dependent on the 
=Fet?,Fet?(Fe??,,.Fet?-Tit")o. content of TiO,. 
@®: Pouillard’s artificial titanomagnetite. 
where the cations inside the bracket 
respresent the 16c position and that outside the 8f position of the spinel structure. 
Hence, the number of Bohr magneton per molecule of any titanomagnetite, 7, will be 


0 10 20 30 TiO; mol% 


Papers gs! n=5(1—x)+4—5(1—*)—4%=4(1—x) (4) 


where the numbers of Bohr magneton of Fe*’, Fet* and Tit‘ are 4, 5 and 0 respect- 
ively. When the content of TiFe.O,, x, increases in the solid solution between TiFe,O, 
and FeO. the linear decrease of saturation magnetization from that of Fe,O; is 


theoretically expected from (4). 
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netite. 


The Curie point may be expected to be proportional to the number of Bohr 
magneton, 2.e. the intensity of saturation magnetization provided that the relation 
established in the case of the magnetism of binary alloy between ferromagnetic and 
non-ferromagnetic elements holds true in the present case of the solid solution bé- 
tween the ferromagnetic and non-ferromagnetic metallic oxides. 

_ As will be seen in Fig. 6, Fig. 7, Fig. 8 and Fig. 9, the linear relations 
between each two among the lattice parameter, the Curie point, the chemical com- 
posion and the intensity of saturation magnetization could be actually obtained as the - 
first approximation for the natural titanomagnetite specimens. If one of the above 
four properties of the specimens of titanomagnetite series is given, we “may be able 
to predict the other three properties with a. fair accuracy. 

(B) ILlmenite-hematite series. 


The intensity of the saturation magnetization of the specimens of ilmenite- 
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hematite series at the room temperature is given in Table II. As will be seen in this 
table, the petrologically so-called ilmenites, No. 1/~No. 4’ specimens, have the small 
value of the saturation magnetizatation such as 0.2e.m.u/gr., being com»arable in 
order with that of hematite,“ while No. 5’ and No. 6/ specimens, which may be 
considered to be situated near the middle of the line DC in Fig. 1, have the fairly 
large value of the intensity of saturation magnetization, 7.e. about 20 e.m.u./gr., being 
comparable in order with that of the specimens of titanomagnetite series. Since the 
source of the feeble ferromagnetic property of hematite may be attributed to the 
parasitic ferromagnetism arised from the local imperfect compensation of the two 
sublattices which will be coupled antiferromagnetically,“ the feeble ferromagnetic 
property of the so-called ilmenites may also be attributed to the same mechanism , 
for that of hematite. However, it seems that the strong ferromagnetic property of 
No. 5’ and No. 6 specimens are too strong to be explained by the same parasitic 
ferromagnetism as hematite. 

The Curie point of this series was also determined from the thermo-magnetic 
curves illustrated in Fig. 4. As will be seen in Table II, these values vary from 
about 100°C to 250°C. No. 1/~No. 4’ specimens have the Curie point around 100°C, 
while No. 5’ and No. 6’ specimens around 250°C. Since the number of the specimens 
treated in the present study is too small, the problem whether or not the Curie point 
of this series changes continuously from about 100°C of ilmenite to 675°C of hema- 
tite will remain for future studies. 

§ 6. Conclusions. . 

From the experimental results obtained in the present study, the intrinsic 
magnetic properties of ferromagnetic minerals contained in igneous rocks were 
clarified to some extent. Especially for the titanomagnetites which were commonly 
found in the igneous rocks, their Curie point, intensity of saturation magnetization 
could be interpreted reasonably -as the magnetic properties of the solid solution 
between TiFe.O, and Fe,O;, a linear relations between each two of the chemical 
composition, lattice parameter, the Curie point and the intensity of saturation magne- 
tization being found as the first approximation. 

Another important point obtained in the present study will be that there exists 
a ferromagnetic phase in ilmenite-hematite series having the peculiar type of their 
thermo-magnetic curve. 

Thus, if the thermomagnetic curve of the ferromagnetic minerals contained in 
the igneous rocks is given, their chemical and crystallographic structure may be 
predicted with a fair accuracy. 

In concluding this paper, the writer wishes to thank Prof. T. Nagata who gave 
him the kind direction and encouragement throughout the study. The writer is 
indebted to Dr. H. Kuno for his useful discussions and suggestions in petrological 
side. His hearty thanks are also due to Messrs. T. Katsusa and J. Osaka who offered 
kindly him the data of the chemical composition of the specimens. 
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On the 27-Day Recurrence of the Cosmic-Ray 
Diurnal Variation* 


By Sekiko YOSHIDA 
(Physical Institute, Nagoya University) 


and Ichiro KONDO 


(Scientific Research Institute) 
Abstract 


27-day recurrence of the cosmic-ray diurnal variation was found 
from the analysis of the data observed at Huancayo during the period 
1936-1940. Similar analysis was done on the daily mean values of 
_cosmic-ray intensity and geomagnetic activities. Among them the 
recurrence of the cosmic-ray diurnal variation was the most persis- 
tent, and that of the daily mean cosmic-ray intensity was the next. 
It was suggested that the most of the cosmic-ray diurnal varia- 
tion may be interpreted as the storm effect. 


1. Introduction 


It is well known [1, 2, 3 and 4] that there are amplitude increase and phase 
advancement of cosmic-ray diurnal variation at the time of magnetic storms or dis- 
tutbances. These facts suggest that not only the worldwide mean intensity, but also 
the anisotropy of cosmic rays are affected by solar stream. Therefore 27-day recur- 
rence of cosmic-ray diurnal variation may be expected. Fortunately, Forbush and 
oan [5] of Carnegie Institution of Washington published the cosmic-ray data from 
Huancayo Observatory of more than 10 years. Harmonic coefficients were calculated 
for every day during this period, to study 27-day recurrence. 

The 27-day recurrence of the worldwide mean intensity of cosmic rays is more 
persistent than that of geomagnetic activity. The persistency of the 27-day recur- 
rence of the cosmic-ray diurnal variation was compared with the persistencies above 


described. 
2. Data used in the analysis 


Analysis of the diurnal variation was done with the bi-hourly cosmic-ray data 
of Huancayo. The cosmic-ray data from Canberra [6] were also used in a part of © 


* Contribution of Geophysical Papers dedicated to Prof. M. Hasegawa on his sixtieth birthday. 
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analysis to obtain the worldwide change of the 
intensity. The geomagnetic data, including the 
daily mean horizontal intensity, the date and time 
4 of the beginning of the magnetic storms and their 
types and ranges, were given from Kakioka Magne- 


tic Observatory [7]. 
In order to simplify the labor of the calcula- 


tion, we made the Fourier analysis of the daily 
cosmic-ray intensity. The procedure of analysis | 


was as follows: 
a) In order to avoid the possible effect of 


long period (>24 hours) variation, 24 hours running ; 
average M, was subtracted from individual bi- 
hourly value x;, thus bi-hourly departure 4; was 


obtained. . 
b) After grouping 4; according to the date 


[aScGn S08 BB 
SOR CSESENGSENSEEE. BeBe 


Littl T RSs 7) 


= ' of U.T., Fourier coefficients were calculated by 
HH ILB.M. Accounting Machines. The 27-day recur- 
HH rence diagram of diurnal amplitude @ is shown in 
Aas es Fig. 1, and some examples of the day to day varia- 
mae tion are shown in Fig. 2. 
WG A 4 1490 
H 
a "| 3. Method of analysis and results 
BEREZ ZAZY 
aay fone? me 
i : it na The study of the 27-day recurrence was done 
7 be : ~ > 
: mee d bi : using Chree’s superposed epoch method. The days 
ial rH of the beginning of the magnetic storms according 
1949 Pb 4 aA to Kakioka report were selected as the epochs (0 
acme . days). The data during the period from 1936 to 


1940 were used to study the recurrence in the sun- 
spot active period. The number of storms used is 
115. The amplitude @ and phase fa, of mean - 
diurnal variation for each storm day was calculated 
from the mean Fourier’ coefficients obtained by the 
Chree’s method. The result are shown in Fig. 3 
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time. 
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Apr. 1937 Jan, 1938 
4 Fig. 2. Two examples of day to day variation of 
I: Daily mean intenstty of cosmic rays 
a: Amplitude of cosmic-ray diurnal variation 


. 
$ To 15 20 25 Jo 


-27 0 +27 +54 day 
Fig. 3 The variation of H, I, a, tmaz, LC.F. and n as the function of storm days. 
For dotted line, see text (§ 3 e) 
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at Kakioka H, LC.F. (International Character Figures of the geomagnetic field), and 
the number of magnetic storms were also summed up by the same method. These 
results are also shown in Fig. 3. 

The following features can be easily seen from Fig. 3. 

a) The amplitude a shows a remarkable 27-day recurrence and its maxima are 
in the neighbourhood of 0, +27 and +54 days. 

b) The phase #,,., Shows a little advancement at the time of magnetic storm 
(+1, +2 days) but shows no recurrence more than their errors. 

c) Daily mean intensity J shows a large decrease after 0 day, and it also 
indicate 27-day recurrence. . . 

d) H, LC.F. and n also show sharp peaks near 0 day, but the persistency of 
the 27-day recurrence are found to be rather small. 

These features are illustrated in Fig. 4 to show their persistency quantitatively. 
This shows the amplitudes of 27-day 
variation of J, H, L.C.F. and a, calcu- 
lated for each 27-day interval having 
centers at —27, 0, +27 and +5! storm 
days, assuming a sinusoidal variation 
for each interval. The amplitudes are 
normalized to the values of 0 day 
interval. The results show the rapid 
decrease of 27-day amplitude at either 
side of 0 day interval in the case of H 
and .C.F. The 27-day amplitude of a, 
however, shows no decrease at the 


“1 0 +1 yt Ragels first cycle, and even at the second 
Fig. 4 The amplitudes of 27-day variations cycle it remains to 0.8 of the original 


Of J.-A; UCL. anda: 


value. Also the 27-day ampli 
The amplitudes are normalized to the : i? coinage’ 
values of 0 day interval. shows slower decrease than that of H 
and LC.F. 


e) The number of storms, », shows little 27-day recurrence, so the 27-day 
variation of @ is not directly connected to the recurrence of the magnetic storms. 
This fact was ascertained by different means. That is also the Chree’s method with 
the same epochs, but the data after the next storms were omitted from the analysis. 
Still, there remains 27-day recurrence as shown by the dotted lines of Fig. 3. From 


these results, we may conclude that the recurrence of the diurnal amplitude is more » 


persistent than that of the magnetic storms. 


f) In Fig. 5, correlation diagrams among Ja, 4J, 4H and IC.F. — shown, 
where 4a, 4 and 4H are the differences between the daily values and the 27-day 
running averages. In these figures, thin lines ‘indicate the values for all days, 


whereas the thick lines indicate the values for non-stormy days, i.e. the days other 
than 0, +1 or +2 storm days. 
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Fig. 5 The correlation diagrams among Aa, AJ, AH and /.C.F. Ordinates represents the mean 
values of Aa for certain bands of abscissa. 

Thin lines: All days Thick lines: Non-stormy days. 


7 The correlation between 4J and da is fairly good, even at the non-stormy 
days, but the correlation with 4H, or LC.F. are very poor at the non-stormy days. 
These results may be interpreted as the result of the difference in the persistencies 
described above. 


g) The results described in § 3 may be summarized as follows. 


The duration or the The persistency of the 
recovery time of a storm 27-day recurrence 
: half decay time half decay time 
“ Cosmic-ray diurnal 
sommes Het longest ~10 days longest ~3 cycles 
Cosmic-ray daily - 
Cotn intennity I longer ~ 5 days longer ~1 cycles 
Geomagnetic activit 
H mee LCF. y short 2~3 days _. short 0.2~0.3 cycles 


4. Discussion 


The cause of the cosmic-ray diurnal variation was tried to be explained by 
several means, for example, Janossy [8], Vallarta and Godart [9] and Dwight [10] tried’ 
to explain it as the result of the existence of the solar permanent. dipole magnetic 
field, but this cannot explain the existence of the diurnal variation at low latitude or 
at lower atmosphere. The possibility of explanation by means of constant solar emis- 
sion of cosmic rays [11], would not be sufficient for the diurnal variation at high 
latitude. 

As these trials for the explanation by means of the parmanent diurnal variation 
seems to fail, we want to propose an idea that there exists no appreciable permanent 
diurnal variation in the primary cosmic radiation. However, when some condition 
around the earth varies, then the cosmic rays are affected and this produces the 
anisotropy or the diurnal variation. 

This idea was drawn from the facts that the amplitude increases very much at 
the time of the magnetic storms [4], and it decays to a very small amount in a rather 
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long time (see for example, the dotted line of Fig. 3). The residual diurnal variation 
can be attributable to the variation of atmospheric temperature (perhaps amounts to 
0.1% amplitudes at low latitudes). 

This hypothesis can be verified if the neutron intensity or the primary intensity 
itself is measured, because they are not affected by the temperature effect. Simpson 
et al. [12] reported that the diurnal variation of neutrons has amplitudes of 0 to 2.5% 
and it varies with time by 27-day period. We have only a small knowledge about the 
existence of the diurnal variation at the top of the atmosphere, but there is no 
evidence about the existence of permanent diurnal variation in the primary radiation. 

The cosmic-ray variation at the time of magnetic storms was interpreted as 
the effect of the cosmic-ray storm [13]. This cosmic-ray storm is characterized by 
simultaneous occurrence of worldwide decrease 4/ and appearance of anisotropy 42 
of the cosmic radiation. The cause of 47 was considered to be the acceleration or 
deceleration of general cosmic rays through some electric field produced by solar 
stream, from the energy dependence of JJ (i.e. latitude, altitude and component 
dependence) [14]. It is also favorable to consider that the enhancement of anisotropy 
was caused by similar electric field due to the solar stream. And the 27-day varia- 
tion of the diurnal variation described above, may be due to the recurrence of the 
cosmic-ray storm. Further, the difference in the persistencies shown in Fig. 4, may 
be interpreted as the time variation of the states of the solar stream. The fresh 
streams can produce magnetic storms and also cosmic-ray storms. However, when 
the stream reaodpear near the earth after one or two solar rotation cycles, it can no 
longer produce large disturbances in the geomagnetic field, but it can produced 
cosmic-ray storms. 
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Anisotropy of Cosmic Rays Observed with 
an Alt-azimuth Telescope 


By Yataro SEKIDO, Sekiko YOSHIDA and Yoshiko KAMIYA 
(Physical Institute, Nagoya University) 


Abstract 


Intensity distribution of cosmic rays over the celestial sphere 
was scanned with an alt-azimuth counter telescope of a narrow 
angular resolution, kept nearly horizontal so that the effective cos- 
mic-ray particles were free from complicated geomagnetic deflection. 
The result suggests the existence of narrow anisotropies, i.e. some- 
thing like point sources. The charge and the average momentum of 
the cosmic-ray particles, with which the narrow anisotropies were. 
distinguished, were observed to be positive and 2 Stérmers respec- 
tively. 


1. Introduction 


Observation of the intensity distribution of cosmic rays over the celestial 
sphere is one of the fundamental subjects in the study on the origin of cosmic rays. 
A counter telescope was designed specially for this purpose. Its angular resolution 
is narrow, so that it may distinguish narrow anisotropies, i.e. something like point 
sources, provided that they exist. The telescope is kept nearly horizontal, so that 
the effective cosmic rays, which pass through a thick layer of the atmosphere, are 
very energitic. Thus, they are free from complicated geomagnetic deflections, so that 
the anisotropy are maintained through the geomagnetic field, and the small geomag- 
netic deflections can be calculated. The telescope is put on an alt-azimuth mount- 
ing, so that a large part of the celestial sphere can be scanned, and the charge and 
the average momentum of the cosmic rays can be measured provided that the 
narrow anisotropies exist. Such a telescope was constructed and the observation has 
been done since December 1951 at Nagaya. ; 


2. Equipment 


The cosmic ray telescope is a 4-fold coincidence set of G-M counters, with an 
alt-azimuth mounting, as shown in Fig. 1. The area of each tray is 868cm?. The 
resolution angle is shown in Fig. 2. The movable ranges are 55°<z<90° and 0°<A 
=<3¢0°, where z is the zenith distance and A is the azimuth. The observation of 
anisotropies was done with a fixed zenith distance z=80°, so as the atmospheric 
depth to be kept constant. In this case the path through the atmosphere was so long 
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== G ' 
Fig. 2 Resolution angle of the telescope 
x: angle between an incident ray 
and the telescope axis 
y: sensitivity. dotted line: geo- 
metrical sensitivity (sensitive 
area). full line: effective sen- 
sitivity, where the zenith 
distance dependence of cosmic 
trays, observed with this 
telescope around z=80°, is 
taken into account. 


Fig. 1 The cosmic-ray telescope. 


that the low energy cosmic rays were filtered off, thus the momenta of the effective 
primary particles being more than about 2 St¢rmers (240GeV). 1588 gr/cm* of iron 
plates were placed between counter trays in order to exclude such low energy parti- 
cles that enter into the telescope after any appreciable deflection in the atmosphere. 
The tray 2 (see Fig. 1) was covered with 5cm Pb to absorb air shower particles, and 
the tray 3 was covered with an anti-counter tray to exclude the effect of air shower. 
The tray 1 and 4 can be slided in their plane. These sliding devices were useful to 
measure the frequency of air showers that could not be excluded by means of the 
ead shield and the anti-counters. The coincidence pulses were registered on a paper 
running 6mm/min. 


q0 3. Preliminary Experiments 


Results of various measurements 
obtained with the telescope kept horizontal 
are shown in Fig. 3. 2f, 3f and 4f are 
frequencies of 2-fold (tray 1-4), 3-fold (tray 
1-2-4) and 4-fold coincidences respectively 
observed with neither the lead shield nor 
the anti-counters. Pb is the frequency of 
4-fold coincidences observed with the lead 
shield only, and Ant is the frequency of 
the 4-fold coincidences observed with the 
lead shield and also the anti-counters. 
H The frequency of 2-fold chance coincide- 


80 


Coincidence per hour 


25 35 4§ Pb Ant 


Fig. 3 Effect due to air showers 
Na: chance coincidences shown in the figure, while it is neglisible 


Ae upaisishoners in the case of 3-fold or 4-fold. The 
HR: nearly horizontal rays 


ncees was estimated, and corrected as 


result shown in the figure indicates that 
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4-fold coincidence, the lead shield and the anti-counters were effective to decrease 
the effect of air showers. The frequency of the air showers that could not be 
excluded was measured after sliding tray 1 and tray 4 out of line (see SL in Fig. 1). 
This frequency was about 1/hour. H in the figure was obtained by subtracting this 
value from Ant. Thus, H is the frequency of nearly horizontal cosmic rays entering 
into the horizontal telescope of a finite aperture. This value H is in accordance with 
that expected from the zenith distance dependence described later. 

At the zenith distance z=80°, where the observation of anisotropy was done, 
the average frequency of the coincidences was about 4/hour. In this case also, the 
frequency of the air showers that could not be excluded was about 1/hour, therefore 
the frequency of the cosmic rays entering along the axis is about 3/hour. 

The zenith distance dependence of the coincidences (air showers subtracted) 
observed with three different thickness of iron absorber, and also the frequency of the 
air showers that could not be excluded are shown in Fig. 4. From this figure it can 
be known that 1588 gr/cm? Fe is sufficiently thick absorber at the zenith distance 
z=80°. The zenith distance dependence are plotted again in Fig. 5, where the curves 
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Fig. 4 Zenith distance dependence 
4a: Observed frequencies, where air 
showers were subtracted. 
4b: Air showers, observed by sliding 
method. —> sec 6 


Fig. 5 Zenith distance dependence 
plots: observation, curves: theory by 
Murayama and others. 
© Normalization point. 
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are the zenith distance dependence expected theoretically by Murayama and others, 1) 
as the solution of a diffusion equation, assuming a source function of z-mesons. Only 
one point in the three curves was normalized with the experiment. Still there is a 
good agreement between the experiment and the theory. 
According to this theory, it is plausible to assume the average momentum of 
the effective primary protons to be about 2 Stgrmers. And, according to this theory, 
it is also plausible to expect the observed flux (3/hour) from the primary protons 
with this average momentum. A direct measurement of the primary momentum is 
described in § 6. 


4. Observation 


The observation was done at Nagoya, where the latitude is 35°N, the longitude 
is 137°E, and the geomagnetic on 


latitude is 25°N. The right ascen- De sch ee ig: 


=. 


sion « and the declination 6 at a 
certain sidereal time, #,;,=0h for 
example, corresponding to a fixed 
zenith distance z=80° and vari- 
able azimuth 0°<A<360°, are 
shown in Fig. 6. The dotted line 
was obtained by neglecting the 


geomagnetic deflection, while the 
full line is the direction of 2 
St¢érmer positive particles at the 
outside of the geomagnetic field 
obtained by calculating their 


i 12 
its. Therefore, the dotted line 
Sates . a Fig. 6 Relation between («, &) and (A, z) at 2=80° 
eeeatep re rer serectasty oF thc and tsig=Oh. dotted line: magnetic deflec- 
2 St¢rmer positive particles com- tion neglected. full line: 2 Stérmer, posi- 
ing from the direction of the full! tively charged particles. 


line. At a certain azimuth A;, a certain declination band (0,) was scanned by obser- 
ving 2 sidereal days (see Table I). Next, the neighbouring declination band (d:) was 
scanned by observing at another azimuth A,=A,+4A (see Table I). Thus, a large 
part of the celestial sphere (—45°<é<+70°, and Oh<v<24h, see Fig. 7) was scanned 
2 or 3 times. Every one point (c, 0) was observed at two different values of azimuth. 
For example, the point P[¢=5h48m é=—5°] was observed at [A=90°, f,:a=0h 00m] 
and at [A=248°, t,:¢=9h44m], as shown in Fig. 7. This combination of the azimuth 
depends on the momentum of cosmic-ray particles. Besides, the relation between 
the solar time of observation f,,, and the right ascension depends on the schedule of 
the observation. The period and the schedule of the observation are shown in Table. 
I, and the relation between /,,, and (a, @) is shown in Fig. 8. 
Every 4-fold coincidence (tray 1-2-3-4) in the series (1), (2) and (3) was regis- 


26 Y. Sexrpo, S. Yosnrpa and Y. Kamiya 


Table I Schedule of the observation 


series Period A n| AA z 


(1) Dec. 1951-Apr. 1952 0°(N)-» 90°(E) -180°(S)}2 | +5° 78.9° +0.1 
(2) Jun. 1952-Sep. 1952 360°(N)—>270°(W)—>180°(S)| 2 | —5° ” 
(3a) Mar. 1953-Apr. 1953 180°(S)—145° 3] =5* 80.0° +0.0 
(3b) Apr. 1953 340°——> 335° A (beret » 
(3c) Jul. 1953-Oct. 1953 280°—» 240° 3.) 6 =5° ” 


12 
Fig. 7 Observed part of the celestial sphere Fig. 8 Positions on the celestial sphere («, 6), 
curve (1): fsig=Oh00m, 0°<A<180° where the observations were done at 
curve (2): tsig=9h 44m, 180°<A<360° tson =20h.f soz =20h is the time of maxi- 
mum intensity in the solar time varia- 
tion (see $5). (1), (2),....: series (1), 
series (2),.... respectively. 


tered. Besides, every 3-fold coincidence (tray 1-2-3 and tray 2-3-4) and also the air 
shower excluded by the anti-counters (tray 1-2-3-4-anti) in the series (3) were regis- 
tered. The average frequencies were as follows :— 

4-fold, series (1)+(2)+-(8) 18842/4691.1 h= 4.02+0.03/h, 

3-fold (2 sets), series (3) 26679/1246.65h = 21.40+0.14/h, 

shower excluded, series (3) 3960/1246.73h =3.18 +.0.05/h, 
As the observations were done continuously during 2 or 3 days at a certain declina- 
tion band, average intensities were calculated for each band. These band averages 
are plotted against the declination, and, on the other hand, 5-day averages are plotted 
against time, as shown in Fig. 9. A part of this variation may be due to atmospheric 
effects. But, the nature is not clear at present. Therefore, no correction was done 
for this variation, but the differences from each band average were subjected to the 
further analysis, thus eliminating the effect of this variation. 
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Fig. 9 Variation with declination (é) and variation with time (ft) 
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(a) —> tow (b) (c) 
Fig. 10 Solar time variation (variation with fs) 10b and 10c: Ist Fig. 11 Two kinds of so- 
harmonics, 10c: 3-fold lar time variation 


@: variation with fsoz, 
4: variation with 0. 


5. Solar time variation and Right ascension distribution 


In Fig. 10a, intensities are plotted against the local solar time ?¢,,,, where ¢,,, is 
the time of the observation, and the direction of the telescope was not taken into 
account. The Ist harmonics are shown in Fig. 10b and 10c. The time of maximum 
intensity is at about 20h in every case. 

In Fig. 11, the Ist harmonics of two kinds of solar time variation are shown. 

- The ist kind (@) is that which was reproduced from Fig. 10b, i.e. the variation with 
Ries rt, the time of observation. The 2nd kind (4) is the variation with 9, the direction 
of the telescope expressed by means of the solar time. As shown in the figure, the 
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results obtained from series (1) and (2) are same as each other for the Ist kind, while 
they are different from each other for the 2nd kind. Therefore, the 1st kind (/,,:) is 
real, and the 2nd kind (0) is an apparent variation caused by the Ist kind. This 
means that the solar time variation is not due to any anisotropy of primary parti- 
cles, but depends on the time of observation. Therefore, it will be attributable to 
an atmospheric effect, but its nature is not clear at present. 


— o 


(a) b c 
Fig. 12 Right ascension pre! Ne, Rane ts with «) = and 12c: Ist 
harmonics, 10c: 3-fold 
X 1 appearent 1st harmonics expectéd as the effect of solar time 
variation. 


In Fig. 12a, intensities are plotted against the right ascension oe. The Ist 
harmonics are shown in Fig. 12b and 12c. The position of maximum intensity is at 
about 12h in the cases shown in Fig. 12b. But, in the case of Fig. 12c, the results 
obtained at various declination are different from each other. However, all the © 
results in Fig. 12b and 12c can simply be interpreted as an apparent variation caused 
by the local solar time variation described above. Expected values are shown with 
x in Fig. 12b and 12c. Therefore the lst harmonics of the right ascension distribu- 
tion must be interpreted as an appearent variation caused by the solar time varia- 
tion. 

If cosmic rays are assumed to be more intense near the galactic plane, then 
2nd harmonic of the right ascension distribution will be expected. The total field of 
the observed celestial sphere was devided into 6 declination zones, and the 2nd 
harmonics of the right ascension distribution were calculated for each of these 6 
zones from the observation of series (1) and (2). Most of the positions of their 
maximum intensities are near the galactic plane as shown in Fig. 13. It is difficult 
to see directly the galactic latitude effect, for the intensity difference from the band 
average has been subjected to the analysis as described in § 4. Therefore, the galac- 
tic latitude effect was roughly obtained by connecting the six 2nd harmonics described 


above. The result is shown in Fig. 14, the intensity near the galactic plane being a 
little higher than that of higher latitudes. 
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Fig. 14 Galactic latitude effect 
(variation with 5), ob- 
Fig. 13 Positions of the maximum intensities tained from the 2nd 
of the 2nd harmonics of right ascen- harmonics of right 


sion distribution. ascension distribution. 


6. Narrow Anisotropy 


The total field of the observed celestial sphere was devided into 1536 small 
squares, each 5° 5°. The intensities at each small square (differences from respec- 
tive band average) were tabulated as the table (1), (2) and (3) for each of the series 
(1), (2) and (3) respectively. For the duration of observation at each square is short 
(from 2.3 hours to 18 hours), it is difficult to find individual point sources or narrow 
anisotropies directly. But the existence of such anisotropies can be discussed through 
statistical arguments. The solar time variation and the galactic latitude effect des- 
cribed above were not corrected in the following analysis. For they are slow varia- 
tions, their effect on the analysis of narrow anisotropy is neglisibly small. 

Assume that any narrow anisotropy does not exist, then the intensity at each 
small square may be considered as a sample drawn from the Poisson distribution 
whose mean value is the respective band average. The squares where their intensities 
are so large that the probability of their appearance are less than 8% were marked 
in each table of series (1), (2) and (3). Let Aq)s@) be the average intensity obtained 
by such calculation that the positions of the squares are determined by the marks on 
the table (1), but the intensities of the squares are read from the table (2). Let Aq) 
be the positions used in the above calculation. Let Ba)» 2) be the average intensity 
obtained by such calculation that the positions of the squares are the surroundings 
(<5°) around Aq), but the intensities of the squares are read from the table (2). 
And, let Ca)5c2) be the averaged intensity at the remaining squares (=10° from Aa) 
read from the table (2). Similar calculations were done also for (2)-(1), (1)>(8), and 
(2)(3). The results are shown in Fig. 15. Intensity at A is larger than those at B 
or C. For the selection of position and the reading of intensity were done on 


Y. Sexo, $. Yosurpa and Y. KamiyA , 


different tables, the statistical 
significance may be checked by 
Thompson’s method. Thus, the | 
significance levelzof the exis- 
tence of narrow anisotropy 
(fixed on the celestial sphere) 
is 0.3%. 

There is little difference 
between the intensity at the 
surroundings (B) and that of the 
background (C), though the 


(1). @)) Sea Se ee el ee angular distance between A 
G) say}average (Pr) GA) Total and B is only 5°, i.e. about the 


Fig. 15 Existence of narrow anisotropies (4-fold) 
A: point in question, B: surroundings 
C: back ground; resolution angle: 5° 


to be not more than about 5°. 


resolution angle of the teles- 
cope. Therefore, the width of 
these anisotropies is considered 


of A, 


Fig. 16 is the result of the similar analysis as that of 


Fig. 15, calculated from the observation of 3-fold coincidence, whose resolution angle 


is about 7.5°. 


In this case, the anisotropy appears to be considerably small. 


This 


result is reasonable, if the width of the anisotropy is not more than 5°. 


In the analysis described above, each pair of the observations made at two 


different azimuth corresponding to one position on the celestial sphere was calculated 
by assuming the geomagnetic deflection of a positively charged particle with the 
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A 

A 
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resolution angle: 7.5° 


x: assumed geomagnetic deflection +10° 
corresponds to positively charged ° 
2 St¢grmer particles - ; 
narrow anisotropies distinguished by - 
assmuing various deflection angles (x). 


Se 


Arisotropy of Cosmic Rays Observed with an Alt-azimuth Telescope 31 


momentum of 2 Stérmers. In this case, the deflection angle is about 10° as shown 
in Fig. 6. Similar analysis as that shown in F ig. 15 were done by assuming various 
values for the deflection angle, including the case of negatively charged particles._ 
The results are shown in Fig. 17, where the notation A or C is the same as in Fig. 
15. This result shows that the anisotropy is missed if the deflection is assumed to be 
too large or too small or negative, but the anisotropy was distinguished if the 
primary particles are assumed to be positively charged particles ‘(probably protons) 
with the momenta of 2 Stérmers (or the energy 240GeV). This means a direct 
mesurement of the charge and momenta of very high energy primary cosmic-ray 
particles. The sign of charge (positive) is in accordance with the general enovulediee 
on cosmic rays, and the momentum (2 Stérmers) is also in accordance with the 
assumed value which was cosistent with the atmospheric absorption as described in 
$ 3. 
4 The number of the small squares marked in the table (1) or table (2) is more 
than a hundred. Now, a new table (1)+(2) was made by adding the intensities at the 
corresponding squares in the table (1) and table (2). In the total field of observation 
(—45°=d=+70°), there were 15 positions, where the intensities in the new table were 
so large that the probability of their appearance by mere chance were less than 1.5%. 
These 15 remarkable positions are shown in Fig. 18, and their distribution with the 
galactic latitude } is shown in Fig. 19. The celestial sphere was devided into 8 zones 
of galactic latitude. The ordinate of 
Fig. 19 is the ratio of the number of 
the remarkable positions to the num- 
ber of observed small squaresin the 
respective zone. Most of these 
remarkable posisions are near the 
galactic plane, though the statistics © 
is very poor. 

Among these 15 positions, there 
were 10 positions, where the intensities 


Fig. 18 Remarkable positions where larger 
intensities were registered in the 
series (1)+(2). 


e : probability of chance appear- pao 

ance S1.5% Fig. 19 Distribution of the remarkable posi- 
©: Pia: S1.0% tions (shown in Fig. 18) with the 
A and B: positions where the galactic latitude 6. dight inside the 
largest intensities were registered figure: number of the remarkable 


in the series (1)+(2)+(3). positions. 
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were so large that the probability of their appearance by nese chance were less than 
1%. Among them, 5 positions are in the part of the celestial sig ice the 
observation series (3) was done. Among the 5 positions, there are 2 agg (A and 
B in Fig. 18), where the intensities in the table (1)+(2)-+(3) were especially esd 
shown in Table II. Assume that there is no narrow anisotropy, then the probability 


Table II 
A B 
h3 
Right ascension « 7h tage 5 ani 
Declination 6 gs a 
Galactic longitude / fh ah 
Galactic latitude D 12°S 


Hours of observation 


Number of coincidence 69 27 


Average intensity in the declination band 43.2/12 pee 
Intensity difference +59% +118% 
Probability w 0.009% 0.014% 
” Wi=Nw 7% 11% 
W.=NwxNw 0.8% 


” 


w to get such a large number of coincidences is very small as shown in Table I. 
However, the number of observed squares in the table (1)+(2)+(3) is N=816. There- 
fore, the probability W., with which this small probability w will happen at any 
two of the squares, is 0.8%, that is to say the significance level with which at least 
one of them (A or B) is considered as a narrow intense position is 0.89. But, the 
probability W;, with which the small probability w will happen at any one of the 
squares, is 7 or 11%, therefore it is not yet significant that both of them (A and B) 
are considered as narrow intense positions. It is desirable that they will be confirmed 
through further observations. 

These two positions are near the galactic plane. Also, the results shown in 
Fig. 14 and Fig. 19 suggest intense positions to be near the galactic plane, though 
their significance could not be proved because of their poor statistics. 


7. Summary and Discussion 


From the observation done with an alt-azimuth counter telescope of a narrow 
angular resolution, kept nearly horizontal, the following results were obtained. 

1) The existence of narrow anisotropies of cosmic rays, fixed on the celestial 
sphere, was shown through a statistical argument. 

2) The charge and the average momentum of the cosmic-say particles, with 
which the narrow anisotropies were distinguished, were observed to be positive and 
2 Stérmers respectively. 

It is well known that the polarization of star light suggests the existence of 
the magnetic clouds, where the strength of the magnetic field is about 107° gauss.) 
The radius of curveture of 2 Stérmer particles in a homogeneous magnetic field of 
10-° gauss is 10cm, i.e.. very small compared with the distance to the nearest star. 
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Therefore it is difficult to interprete these narrow anisotropies as the point sources 
of cosmic rays, provided the solar system is in a magnetic cloud. 

Two of the authors”) have suggested Crab Nebula as a point source, but no 
anisotropy was observed by Barrett and others”. The statistical accuracy of the 
result obtained until now by this alt-azimuth telescope is insufficient to discuss on 
this discrepancy. 
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An Analysis of Electric Field after Lightning Discharges* 


By Yuichi TAMURA 
(Geophysical Institute, Kyoto University) 


Abstract 


Analysis is based on a following set of simple assumptions. 
A thundercloud is represented by a vertical dipole of positive polarity. 
Electrical conductivity of the atmosphere increases exponentially with 
height, and not altered by presence of thundercloud. Electromotive 
force in thundercloud is represented by a constant current generator. 
It is taken that an electric field consists of three parts. The first is 
field due to charge in each pole which increases at a rate depending 
on electrical conductivity of the air at a height of pole; the second, 
due to space charge which grows with charges in poles; the third, 
due to free decaying space charge which is released after a lightning 
discharge. Characteristic dissimilarity of the form of field recovery 
curve, observed after a near discharge to that of distant discharge, © 
is well interpreted by this analysis. a 


1. Introduction 


One of the most interesting features of the records of electric field during 
a thunderstorm is the “recovery curve” of the field which follows a lightning dis- 
charge. Among many types of the field recovery curve (see references [1]2]3]4)), 
the simplest approximates to exponential form which is usually seen when a discharge 
is sufficiently distant. The characteristic feature of this simple type has been 
regarded as an important clue yielding an information about the manner of separa- 
tion of charges in thunderclouds. Suggestion relating to such simple recovery curve 
has initially been put forward by Wilson [1] in the following manner. In a thunder- 
cloud there presumably exist carriers of opposite electricity having sizes distinctive 
to the sign of electricity they carry. In the cloud, being originary neutral, vertical 
separation of charges will begin owing to the relative displacement of these carriers 
under gravity. In such circumstances the rate of accumulation of charge is not 
constant, but becomes less as the accumulation progresses, because two kinds of 
retardation may act; one as an effect of electric field in the cloud against further 
accumulation of charges, and the other as dissipation of charges by conduction. 
Assuming the rate of accumulation under no opposition be constant and each retard- 


* Contribution to Geophysical Papers dedicated to Prof. M. Hasegawa on his sixtieth birthday. 
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ing effect be proportional to charges accumulated, it will be seen that the charges 
have to increase exponentially with time until a lightning discharge occurs. Further, 
if the distance between oppositely charged portions remains unchanged, the charges 
themselves will be proportional to the electric moment (defined as 2 x charge x distance 
between opposite charges) of the cloud. Since an electric field intensity, fine weather 
field being put out of consideration. observed on the ground at sufficiently remote 
distance from a storm, is proportional to the electric moment of the cloud, the form 
of recovery curve of the field is also exponential having a time constant equals to 
that of charge accumulation. Wilson’s suggestion has further been emphasized by 
Wormell [2]. As far as the elementary process of charge generation in thunderclouds 
is concerned, theory of Simpson differs from Wilson’s. They, however, are in agree- 
ment with respect to the manner of charge separation [5]. 

Meanwhile, there is a point of view that the atmospheric electric field be con- 
sidered as field accompanying electric current rather than a static field. This is 
substantial if the spatial variation of electrical conductivity of the atmosphere is 
taken into acount. In recent years Kasemir [6] has discussed various atmospheric 
electrical phenomena in such atmosphere and has given a remark relating to the 
recovery curve of field after a lightning discharge ; Holzer and Saxon [7] have estimated 
the electrical conduction current in the vicinity of thunderstorm under equilibrium 
condition, assuming the atmospheric conductivity increases with height. 

In the present paper it is attempted to analyse the electric field due to a 
thundercloud after a lightning discharge under assumptions that a simple dipole cloud 
is imbedded in the atmosphere conductivity of which, same with that of Holzer and 
Saxon, increases with height, and to go one step further in an interpretation of 


recovery curve of field. 
2. Characteristic Features of Recovery Curves 


Records of electric field at the ground during a thunderstorm of positive 
polarity in which, as often experienced in the later stage of storm activity, intracloud 
lightning discharges of similar scale occur repeatedly is schematically represented in 
Fig. 1. At time f, é&...., the #+1th, 
n+2th,.... discharge occured respec- 


tively. The upper figure corresponds Near storm 
to near storm and the lower figure % wu 
corresponds to distant storm. In the Be ae 


former case discharges occur when the Distant storm 
electric field is in a diminishing stage bee ee ee J aids Tete ye 

udden field changes accompanyin 
snipe ; : 2 re hes nti Nt2 = 1H3 NA 
discharges, are positive; while in the : : 5 = 

i e time : 2 3 4 

ethan, ree cliacsarge in ig ei, Fig. 1 Schematical representation of field 
after the field has attained its station- iecovery cirve. 
ary value and sudden changes are 
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Fig. 2 Examples of electric field of thunderstorms (observed at Kyoto). 
The upper record is for near storm (16"10™-16"34™ Aug. 16, 1944); figures in the 
record are distances (km) of flashes measured from lightning-thunder interval. The middle 
is for intermediate distant storm (1951™-2003™ July 25, 1950); thunder was audible. 
The lower is for distant storm (20"36™-21500™ Sept. 1, 1943); no thunder was heared. 


negative. The definitions of the near the distant storms are somewhat ambiguous, 
it may rather be distinguished by the forms of recovery curves. However, accord- 
ing to the experiences of the writer, with allowance of certain extent, the discharges 
within 7km and beyond 15km will be called the near and the distant storm res- 
pectively.. At intermediate ,distances from discharge there must be a border line 
passing over which a sudden field change reverses its sign. The field recovery 
curves at these distances are somewhat complicated. In Fig. 2 examples of record 
of near, intermediate distant and distant storm obtained by means of the rotating 
collector deviced by Hasegawa[8] are reproduced. Although they are not the simul- 
taneous records of same flashes, the characteristics of recovery curves seenare most 
familiar ones to each distance. Therefore the correspondence shown in Fig. 1 is 
conclusive. 

It is of interest to examine, by employing a simple model of thundercloud and 
the atmosphere of uniform electrical conductively, the results of observation men- 
tioned above may be interpreted or not. One of the most simple model of cloud is 
vertical dipole of positive polarity, opposite charges being equal in magnitude and 
they are distributed uniformly, for example, in spherical volumes. In this case 
electric field intensity on the ground at sufficiently distant from a storm, is pro- 
portional to the dipole moment. But, when the polar distance is assumed to be 
unchanged, electric field intensity is everywhere proportional to the charge itself. If 
after a lightning discharge, the charges accumulate exponentially with time, corres- 
poalingly after a sudden field chang> which is positive ani negative for near 
and distant discharge respectively, the field will recover showing exponential 
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form with a time constant equals to that of charge accumulation. Since these cir- 
cumstances do not agree with observation in which recovery curve for near storm 
-is clearly dissimilar to that for distant storm, the simple picture employed above 
must be abandoned. However, one of the factors affecting the recovery of the field 
will be mentioned here. It is quite probable that below the base of cloud, positive 
space charge may appear by point discharge from the earth’s surface where a strong 
negative field prevails. But this space charge may influence, as it is disappointed, on 
recovery curves to deviate much more from observed ones than when there is no 
space charge of this kind. 

On the other hand, the most striking feature of field recovery is that of distant 
‘storm, that is, as already mentioned, a discharge occurs some time after the field at 
the ground has attained its stationary state. In case of the electric field within 
cloud grows directly proportional to charge in poles, as being assumed here, discharge 
will occur when the field intensity reaches the critical value, or in other word charge 
in poles accurhulates to the limiting amount. It is quite likely that this condition, 
in general, be attained when accumulation of charges is going on and not at the state 
when accumulation of charges have already ceased. Conditions of initiation of 
lightning discharge may be quite delicate, that is, heterogeneous distribution of charge 
may strengthen a local electric field sufficiently large to initiate a discharge. even 
when the general tendency of the field is still unfavorable for a discharge. However, 
it is very difficult to expect that discharge occurs some time after the electric field 
intensity, or the accumulation of charges has reached the stationary state. For, from 
‘statistical point of view, general circumstance of the field will determine the chance 
of initiation of discharge allowing a certain fluctuation depending on the heterogenuity 
in the field. Therefore, as far as atmospheric conductivity is assumed uniform, a 
serious difficulty lies in employing the simple model of cloud. If, however, at the later 
stage of charge accumulation, an assumption is taken that the polar distance varies 
so as to keep the electric moment be constant instead of the polar distance be cons- 
tant, then electric field intensity at remote distances from a storm will have attained 
its stationary value in spite of the field within cloud is still growing towards the 
sparking limit. Since the characteristic feature of field recovery curve for. distant 
storm may be interpreted to some extent. But this tentative way of thinking is not 


developed here. 
Meanwhile, in many years past, an important factor relating to escape of 


charge from thunderclouds was pointed out by Wilson [1]. He stated that in a di- 
polar cloud, rate of dissipation of upper charge will be larger than that of lower 
charge on account of the atmospheric conductivity is larger at the top of cloud than 
at the bottom; especially he suggested that the strong electric field at the top of 
thundercloud (being assumed to be positive polarity) may drag down negative ions 
from the upper conducting layer to make the conductivity of the atmcsphere much 
greater than fine weather value. Although recent observations carried out by Gish 
and Wait [9] showed Wilson’s suggestion does not hold and that atmospheric conduc- 
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tivity over the top as well as in the vicinity of thunderclouds do not differ from that 
of fine weather, it may certainly be said thundercloud is immersed in the atmosphere 
whose conductivity monotonously increases up to a layer higher than the top of 
thunderclouds. Also, Wormell [3], in his detailed discussions on the recovery of the 
electric field after a lightning discharge, has shortly remarked that dissipation of the 
charges will not in general be the same for both charges. But he did not further 
discuss its effect on the recovery of the field. ’ 

According to the writer’s opinion, one way possible to interpret the charac- 
teristic features of field recovery, which have been emphasized above, is to take 
account of nonuniform atmospheric conductivity instead of unform, as Gunn [10] as 
well as Holzer and Saxon [7] employed in their studies on electrical phencmena of 
thunderstorm. In these circumstances, space charge outside the poles should play 
a role on the electric field. In the following paragraphs an analysis of the electric 
field after lightning discharges is attempted from this stand point. 


3. Model of Thundercloud and Assumptions 


Analysis is based on a following 
simple model of thundercloud and fun- 
damental assumptions. In Fig. 3, A and 
B represent positive and negative poles 
of a thundercloud respectively forming 
vertical dipole fixed in the atmosphere. 
A' and B’ are the image of A and B 
respectively. Separation of electricity is 
. assumed to be constant rate, that is, 
electromotive force is equivalent to a 
constant current generator terminals of 
which are A and B. Electrical conduc- 
tivity of the atmosphere is assumed 


. Fig. 3 The model of thundercloud and 
to increase exponentially with height coordinate system 


and not altered by the presence of thundercloud. 


Coordinate system employed is cylindrical (7, z); upward directing axis of the 
dipole is taken as z-axis and its intersection with ground is taken as the origin, 
centres of A and B being (0, z,) and (0, 2.) respectively. apa i of A and B are 
assumed small relative to 2, and 2. 

Effect of the upper conducting layer which is found as high as several tens 


kilometers above ground is ignored, and the normal fine weather electric field is also 
ignored. 


4. Analytical treatments 


As the equation to be solved is, as shown after, linear form it may be treated, 
conveniently, in the case of single pole; when two poles are present the solution is 


| 
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given by superposition of solutions for each pole. 

At first, total charge Q, the function of time f, in the positive pole will be 
considered. Let J).be total current supplied to the pole, and being assumed constant ; 
total conduction. current from the pole to the atmosphere will depend on electric 
field intensity as well as on conductivity at the surface of the pole. But, 
as the dimension of the pole is assumed small relative to its height, the effective 
conductivity may be represented by o; which is the conductivity at the height 2. 
Then, putting 4:=470, and Q,=< 


Q=Qie-*1* + Q,(1—e71*) (1) 
where Q; and Q, are initial and stationary value of Q respectively. 


Electric field intensity H,, due to Q, at any point P(z,7) in the atmosphere 
excluding the region of the pole is represented by its components 


2, 2t 
Pain Ep .= Q,( e252) 


: (2) 
| ae =O)9 (e- a) 
where R={(z—z2)'+7°} and R’={(zt+2,)2+77}12 
Similarly, field intensity E,;, due to Q; is 
a= Oi( Set) 
: (3) 
Exi,,= =@Q; ra ae 3) : 
Therefore, field intensity #, due to Q is 
E,=Eqneo'+E,(1-e"1) (4) 


Let i be conduction current density, # electric field intensity, o space charge 
density and o conductivity of the air. Then following relaticn hold at any point 
P(z,r) in the atmosphere except the pole under consideration 

divi=— ue (5) 
Since t=cF ard divE =4n7p 
} dio( + iE) =0 (6) 
where 4na is replaced by 2. 
Special integral of equation (6) is 


OE 7 

OE ‘+ AE =0 (7) 
Equation (7) is satisfied by 

E=E,.e** > CBi 


showing that initially existed field decays with en, 
General integral of equation (7) is 
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CE + 1B = 1B, r, t) (9) 
where EF, is arbitrary vector satisfying the relation div iE.,.=0, and must ke deter- 
mined so as to F satisfies the initial and final ccnditicns as well as boundary con- 
ditions in the problem. Ce conditions are as follows, 

(i) at <=0, H=0 and op ah Bow 
this conditicn means if there is no initial field, consequently no charge be 
found anywhere, the initial rate of field change must be determined by the 
initial rate of charge accumulation in the pole. 

(ii) at tow, EE, (FE, is the final value of #) 

(iii) at any time, H must be vertical at the ground and vanishes at infinity, 
and the relation Q=Q,(1—e->1') always hold at the pole. 

Then EF, is found as 


AE =iE,l—e 1") +4, F267" (10) 


where £,, is given by equation (2). 
Frem (9) and (10) 


SEB =1E-e"') +h Be ion wef a 

Equation (11) is satisfied by 
E = EK, [des ty Atent —e"™)! - Ey "5 (e7>1¢ —e->*) (12) 

= 


If at ‘=0, H=£;, solution of equation (11) is 


E=E,e™“+E fq gaye AU gene e-méyl E,. Ai (e->1f— e->*) (13) 
: a is J Aaa 
putting H,=E,,+E,., where E,, is the field intensity due to charge in the pole at 
staticnary state, and F’,, is the field intensity due to space charge in its stationary 
distributicn ; equation (13) is conveniently expressed in the other form as 


E=E,,\—e-'))+ Bp 1-eN' +4 A (ete e™)| 4 Bee (14) 
1 “ 


Equation (13) or (14) is general result of the present analysis indicating, as equation 
(14) shows, that the electric field consists of three parts; the first is field due to 
charge in the pole, the second is field due to space charge, and the third is free 
decaying field. EK, the final value of the field intensity is obtained from the relation 


div AE.,=0 (15) 


Using the fundamental assumption on the conductivity of the atmosphere, that is 
ao=o,e"**, equation (15) becomes 


dw E',+2kE,,,=0 (16) 


where £.,. is vertical component of E,. 
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This is equivalent to the equation which Holzer and Saxon [7] have derived and 
solved, namely 
: | 
prot+2k oe =0 (17) 
where ¢ is potential of E,. 
The solution of this equation fit to the present aim is 


KR -kR’/ 
iin ) (18) 


then the electric field intensity E, is easily obtained. 
Putting z=0 and A=Ap, (Ao =4 770») in (13), electric field intensity at the ground 


is given by 
E=E;- t+ BS e-*1 hes a =e te “| Exs7 Ai =e (eit eee) 159) 
with ! 
BREN a naoonst te zt (20) 
js 7 Be (21) 


where R=(r?+2z,°)'/? 
See coponcins to the equation (14), equation a becomes 


E=E,e,0-—e1 t)4 Eos) 1—e- ss ee say Nis e~*o )) HE, re (22) 
{ 
5. Numerical Calculations 


Numerical values are taken as: 
2,=7km, z,=5km the heights of positive and negative pole respectively ; 4; =1.8x 10-3 
esu, ¢>=1.2 x 10-* esu, the conductivities of the atmosphere at the height 7km and 5km 
respectively, assuming ¢,=410-'esu and k=0.11 per km: and: Q, are 80 coulombs 
and —120 coulombs for upper and lower pole respectively, assuming J,=1.81 
ampere. 

Calculations are confined to the electric field at the ground, because no obser- 
vation is to be compared with except that. 

(a) Stationary State—Electric field intensity at stationary state F,, which must 
be attained if no lightning discharge occurs, is calculated from equation (20); also 
E,,, the field intensity due to charge Q,, is calculated from equation (21); then E,, 
=E,—E,, is the field due to space charge in stationary state. These magnitudes at 
various distances from the origin are given in Table 1. 

It will be worth while to note that space charge affects to strengthen E,, at 
near distance and weakens it at remote distance, especially this effect is very large 
at great distance. 

(b) Rising state—Electric field intensity in rising state F is calculated from 
equation (22). Examples of these magnitudes varying with lapse of time are given 
in Table 2, and in these calculations initial field is assumed absent. The correspond- 


ing magnitudes of FE, are also given in the Table. 
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Table 1. Stationary electric field intensities at various distances. 


Q,=—120 coulombs at z22=5km 


Q;=80 coulombs at z3=7km 
ry (km) Es (v/m) Es (v/m) Eqs (v/m) E, (v/m) 
oF Ri a? “204 < 10? 520 x 10 —864 x 10? —13410* 
4 192 « 10? 32410? —411 x10? —601 10" 
8 838 x 10! 122 x 10? —129x10?: —161X 10° 
12 376 x 10! 444 x10! —492x10! —496 x 10! 
16 189 x 10! 17410! —229x10! —179x10! 
20 107 x 10! 744 —124x10! —721 
25 576 273 —652 — 260 
30 344 110 —384 —102 
40 150 204 x 10-1 —165 — 185% 10-* 
50 783 x 10-1 476 x 10-* —850 x 107! —384x10-* 


Table 2. Electric field intersities (volt/meter) versus time (second) 


| at r=0km | at r=50 km 

3a a ee ee BE Ey, | 2, | &- | & pe 
0 0 0 0 0 0| 0 0 0 0 0 
10| 5970] 6130 |~12100 |~12400 |— 6300} 15.90 | 15.38 |—11.90 |-11.41| 3.97 
20} 10700} 11200 |~22600 |—23300 |-12100 | 28.52 | 26.83 |—22.16 |—20.94| 5.89 
30 | 14500 | 15400 |—-31500 |-32900 |-17500 | 38.62 | 35.60 |—31.03 |—28.67| 6.93 
40 | 17500 | 19100 |~39100 |—41600 |—22500 | 46.67) 41.59 |-38.51 |-34.29 | 7.30 
50 | 19900 | 22100 |—45800 |—49500 |~27400 | 53.09! 45.81 |—45.05 |-38.80 | 7.01 
60 | 21800 | 24800 |—51400 |—56300 |-31500 | 58.18 | 48.47 |—50.58 |—42.22 | 6.25 
70 | 23400 | 27100 |~56400 |~62600 |~35500 | 62.25} 50.18 |~55.42 |-44.86 | 5.32 
80 | 24600 | 29100 |—6c600 |—68100 |-39000 | 65.46 | 50.96 |-59.59 |-46.78| 4.18 
90 | 25600 | 30800 |—64100 |~73000 |—-42200 | 68.12 | 51.21 |-63.16 |—47.96 | 3.25 
100 | 26400 | 32400 |—67200 |—77500 |—45100 | 70.16 | 50.74 |-66.22 |—48.70| 2.04 
110 | 27000 | 33700 |~-70000 |~81700 |—48000 | 71.80 | 50.00 |-68.85 |—48.85 | 1.15 
120 | 27400 | 34800 |~72200 |—85300 |~50500 | 73.13 | 49.08 |—71.06 |—48.76 | 0.32 


Ej+ represents the field due to the positive pole, Ey represents the field due to the positive 
pole with space charge accompanying it; E,. and E_ are the corresponding fields for 
negative pole; E(=E,+£_) is the resultant field. 


The following will be noted. If the case of single pole is considered, the effect 


of space charge to the electric field is not neglegible even in the earlier stage of. 


charge separation; while the case of two poles forming a dipole, space charge does 
not affect much to the electric field. Therefore, it may be. said, at least in the 


earlier stage of charge separation, the electric field due to a dipole is determined as: 


if the dipole were in the atmosphere of uniform conductivity. It is important, however, 


to remember that each pole has charge of opposite sign, but in general of different. 


magnitude. 


(c) Field recovery curves—To obtain field recovery curves which are the 


graphical relation of electric field intensity versus time after a lightning discharge, , 


numerical values given in Table 2 are available. 


However, it is neccessary to make 
some assumptions as follows :— ; 
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(i) lightning discharges occur between charges in the upper and the lower 
pole; no cloud-ground discharge nor lightning towards atmosphere occurs, 
(ii) by a lightning discharge half of total positive charge in the upper pole 
and corresponding negative charge in the lower pole are neutralized, 
(iii) space charge outside the poles does not concern a lightning discharge, 
(iv) lightning discharges occur with a constant time interval, this is the time 
taken for accumulation of positive charge to a specified quantity. 
Together with above assumptions, following remarks become important. It can 
be said in equation (22), the electric field intensity due to space charge, namely 
B= Eyles Hty4Ge 
in pole, namely E,=£,,(1—e-*:'). Therefore, by a discharge occured at /=r, if 
E,=E,,(1—e-*1') changes suddenly to E,/=E,,(1—e-*1*0), where r>1)=0, the regene- 


e~*1'— e~*o')', always have to accompany that of due to charge 


ration of FE, thereafter will be represented by 
E,= E,/e7*1¢-*) + EF, (1— e741 -*?) = FE, ,(1— e781 0#*#-*) (23) 


Consequently, at a instant of the flash, E, is divided into two parts; the one is 
accompanying FE, and the other is free decaying part. They are 


Bia; /1— e710 ilies pe Ai ee ~T4t)__ Q-AglT orae0y| (24) 
rae 
YF Fo in emia (25) 
where Blok. fee L Ai (e741* —@-A17 0 — e- rot + erro)! 
? | A1—Ao 


putting ¢’/=i—r, after the discharge E is represented by 


E=E,,i—e1 Co") + Es} 1—2-4:1C ot *? 4. oy (e710 ot?) —@ erature} 
LE et Ble Rot! ‘ (26) 


In practical calculations it is convenient to select the instant of a flash as a new 
origin of time, consequently initial conditions must be selected such as, at ¢/=0, 
E=E,e>0° +E. 

Now recovery curves of field will be found as illustrated in Fig. 4 (this figure 
corresponds to the case of the field at r=50km, the lightning interval is about 74 
seconds and there is no field initially). Curves £,, and E, represent the electric 
fields due to the positive pole and that due to the positive pole with space charge 
accompanying it respectively (see equation (22)); curves E,- and E_ represent corres- 
ponding fields for the negative pole; curve E represents the resultant field (E,+£_). 

It is assumed a half of positive charge in the upper pole and corresponding 


negative charge in the lower pole are neutralized at t=100sec. Then, E,, changes 
from a to c, correspondingly E,. changes from a toy T= 5x83 xac, referring 
Table 1). Consequently E changes from e to f (¢f=ac—ay). After the discharge each 


field which. was formerly represented by EF;,, E,, E- and E_ will grow along £’),, 
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40 80 120 10 200 240 
Fig. 4 Method of drawing recovery curves. 


E’,, E',- and E’_ respectively in the same manner as if the formers start from m, 
n, € and y respectively. The resultant field E’ (=E’,+£’_) will grow starting from 
g. The magnitude of electric fields released by the discharge of positive and negative 
charges in the poles are cd—ab and of—yé respectively, therefore total released field 
is cd—ab+oB—76 which equals to —gf. Then the magnitude —gf is the initial value 
of the free decaying field released by the first discharge. The resultant of E’ and 
the free decaying field is the recovery curve after the first discharge as represented 
by £ (after 100 sec). 

The second discharge occurs at =174sec. At this instant, similarly as in the 
case of the first discharge, E’,, and E’,. change from a’ toc’ and a’ to 7’ respectively. 
Consequently E’ changes from e’ to f or E changes from 7 to j’ (e’f=7j'). The 
curves E”,,, E”,, E”;. and FE”. shown after the second discharge have the same 
meaning with E’,,, E’,, Eg. and E’_ respectively. The resultant E”’ (=E”’,+E”_) 
grows starting from g’. Therefore the magnitude of —g’j’ is the initial value (at 
t=174sec) of the free decaying field which consist of two parts, the one is released 
by the second discharge and the other is remainder released by the first discharge. 
The resultant field of E” and the free decaying field are represented by E (after 
174 sec) showing the field recovery curve after the second discharge. 

By similar procedures the field recovery curves after the subsequent discharges 
will be obtained. . 

Using this graphical method, recovery curves for various values of y are 
obtained and shown in Fig. 5. In this case it is assumed at #=0, electric field was 
absent and at f=90sec. the first discharge occured neutralizing 35 coulombs (=% x80 
x (1—e7#2x1-8x107 °x00) coulombs). Thereafter the discharges repeated with the constant 
interval of 67 seconds. Curves for =0 and r=20km are representatives of near and 
distant storms respectively. Curves for y=8 and 12km are examples of storm at 
intermediate distances. These curves well represent the characteristic features of 


rss 
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second 90 157 224 291 
Fig. 5 Field recovery curves at various distances (calculated) 


those in Fig. 1 as well as Fig. 2. 
6. Discussions 


Relating to the mathematical analysis and numerical calculations performed in 
the previous paragraphs, the following remarks will be necessary. 

(i) The recovery curves in Fig. 5 were obtained assuming half of total posi- 
tive charge (and corresponding negative charge) be neutralized. In the case of total 
positive charge are neutralized, important features of the curves will not alter. Also, 
it is obvious the form of the curves does not depend on the magnitude of supply 
current J,. On the other hand, the characteristic features will depend on the interval 
of lightnings. For more frequent discharges than asstimed in the example, charac- 
teristic features of field recovery for a distant storm becomes less marked, namely, 
discharge will occur when the field is in rising state. This matter is also the same 
in actual storms having very frequent discharges. 

(ii) It has been assumed the supply current J, is constant. If this current 
diminishes owing to the retarding effect of the electric field within cloud, instead of 
equation (1), Q will be expressed by 


Q=Q,e*1' +Q,(1—e-O1**) (27) 


where Q,= F fe and y is a positive quantity depending on the retardation. This 
modification, however, does not affect the course of the analysis remembering 4, be 
replaced by 4:+7 except for free decaying term. For example, equation (22) is 
replaced by ; 


Ay+r g rei e 
E=E, (1-6-0!) + Ey,} 10-010 Liew ely Ott 6 rot) + Bye yt (28) 


where E,, and £,, are , 2 times those in equation (22). Resulting recovery curves 
1 
will have similar features with those shown in Fig. 5. 
(iii) In case of a simple dipole cloud is imbedded in the atmosphere of uni- 
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form conductivity, as previously Wilson [1] pointed out, following relation will hold 


(37) 

ot t=0 3 Ja (29) 
E Q 

is initial rate of field recovery, Eis a sudden field change accompa- 


Oe 
where Ge : 

Je * 
nies lightning, J) is the supply current in a cloud and @ is the charge destroyed by 


lightning. Since the left hand side of equation (29) is measurable, as has been done by 
Wilson [1], Schonland [2] as well as Wormell [3], Jo/@ will be estimated by the obser- 
vations of electric field on the ground. 
According to the present analysis, corresponding relation is 
(Se), tr 40k I, 
fee ___ = (30) 


PRIS ~@Q 


Equation (29) as well as equation (30) are obtained assuming initial accumula- 
tion of charge is absent. If initial charges are present these relations will become 


more complicated form. 
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Rock Magnetism as a Temperature Indicator 
By H. MANLEY 


Abstract 


The magnetic properties of a suite of Scotch dykes are given 
briefly and the results afford a good correlation with their degree 
of metamorphism. The intensity of remanent magnetisation falls 
rapidly with increase of metamorphism. The unaltered dyke material 
described here is believed to be a case of remanent magnetism in 

é the oldest rocks yet mentioned in literature. 


Introduction 


This work describes results obtained for a set of metamorphic dykes in 
Scotland and in which an attempt is made to correlate the rock’s magnetic properties 
and their metamorphic history. They are of interest in so far as these rocks seem 
to be quite the oldest known ones to be mentioned in the literatures as proved to 
possess appreciable remanent magnetisation. 

Cubes were cut from samples of a suite of Pre-Cambrian Pre-Torridonian 
dolerite dykes of the west coast of Sutherland and Ross, Scotland. The four dykes 
considered here comprise a series grading from unaltered injecta to strongly regional 
heat-metamorphosed hornblende-schists and amphibolites. The geological implication 
is that the metamorphism was largely the result of strong heating with Plat ialy 
modest pressures rather than the converse and this is of some significance in the 
interpretation of the results. (1). 

It was intended to take susceptibility and remanent magnetism measurements 
prior to making magnetic traverses across the dykes and collecting oriented samples, 
by means of which the direction of magnetisation of the dykes could be determined. 
Brief petrological and geology notes are given below, on the dyke and dyke rocks. 
The geological data is due to Drs. J. Sutton and J. Watson of the Geology Dept., 
Imperial College, and the petrology regarding the magnetic minerals was confirmed 
by the author under their supervision. 


Petrology and Geology 


Dyke S-445. is an unaltered dolerite dyke which cuts pyroxene gneisses. 
Maenetite is present in this rock up to, say, 8%, and in very large crystals, quite 
the largest for any of the magnetite in the rocks studied in this research. 

Dyke S-2021. is a metamorphosed dyke containing 2% of sphene and rare 
ilmenite. Ta2 magnetite of th> original rock has been altered by the regional | 
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metamorphism to sphene and hornblende. 
Dyke S-243. is a hornblende schist derived from the original dolerite which 
has been strongly heat-metamorphosed.. It contains some sphene and very little 


magnetite. 
Dyke S-100. is an amphibolite derived from dolerite which has been magnetite. 


It may be assumed that these dykes were originally similar dolerites when they 
were injected in Pre-Cambrian times. Since then, however, and the geological im- 
plication is that this occurred relatively soon after the injection, the three latter sets 
of dykes suffered a regional metamorphism which has completely altered their 
texture and mineral content turning them into hornblende-shists and amphibolites, 


the latter being the final stage in the process. 
The unaltered dyke outcrops on the edge of Coch Torridon, Ross-shire, (see 


Fig. 1), and the other dykes are from Loch Inchard and Loch Laxford in Sutherland, 


close by. 
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Magnetic Resuits 


TaBLe 1 
As Cut 
Cube Kn. R alg 
S-445a 7,600 2,100 0.51 
1D 7,500 1,900 0.48 
S-2024a, b 65 20 0.6 
S-100a, b 60 20 0.6 
$-243a, b 80 10 0.2 


Remarks on ‘ As Cut’ Results 


Two samples of each rock type were tested, the results for the altered dykes 
being so similar that the mean values are given. The unaltered dyke, S-445 has a 
very high susceptibility and strong remanent magnetism, whereas the metamorphosed ° 
dykes have low susceptibilities and almost non-detectable remanent magnetism. 

The unaltered dyke rock has a susceptibility at least twice as great as that for 
any other rock tested in this research, and this is consistent with its high magnetite 
content. The remanent magnetism is moderately large, considering the extreme age 
of this rock which is at least 600 million years on Holmes radio-active scale. It is 
interesting to compare the mean T value of 0.5 with that of the Whin Sill, which is 
2.5 a similar type of English rock and only about 1/4 as old. This is the first time 
that any rock older than Paleozoic has been found to possess remanent magnetism ; 
indeed Koenigsberger stated in 1938 (2) that no rocks in Central Europe older than 
Mesozoic possessed remanent magnetism. One must conclude that the terrestrial 
magnetic field existed in Pre-Torridonian times and it would appear, before consider- 
ing the results of the heat-treatment, that it would have been of the. same oder of 
magnitude as at present. 

The three metamorphosed dykes present remarkably similar magnetic proper- 
ties. Their susceptibility is low but readily detectable, while their remanent magne- 
tism is at the lower limit of sensitivity of the apparatus. It would seem that the 
metamorphism to which the dykes were subjected has destroyed all but 1 or 2% of 
the magnetite—a result which is in excellent agreement with the petrology. This 
would seem to be, also, some evidence that rocks found today with strong magnetic 
properties have not been subjected to such regional metamorphism. This is of 
particular interest since there exists a considerable body of literature on the magnetic 
properties of igneous rocks (3 & 6), and most of the rocks studied have been basic intru- 
sions and injections of basalts and doleritic—andesites, similar to this rock type. 
The evidence and assumptions involved in the major interpretation of the magnetic 
properties of ancient rocks such as these dykes is considered in the final summary, 

It must be presumed, tentatively, that these results are strong evidence for 
the existence of the terrestrial field 600-800 milion years ago. 

The results obtained after heat-treatment are :— 
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TABLE 2 


After Heating. 


Cube Kt ».4 EP Ss Q Z 
S-445a. 6,050 0.80 7,900 3.8 2.4 4.75 

Rees 5,650 0.75 6,100 32 2.0 4.2 
S-2024b. 65 1.0 40 2.0 1-2 2.0 
S-100a, b 60 A) 35 Le 1.0 ules 
$-243a. 125 1.6 190 19 2.8 12 

Ds 110 1.4 135 14 23 10 


Remarks on heated Results 

The unaltered dyke rock, S-445, has lost 20% of its susceptibility, presumably 
due to decomposition of its magnetite content, and gained remanent magnetism by 
a mean factor of 3.5 giving a resultant average Q of 2.2. It is interesting to compare 
the effects of this heat-treatment with the natural metamorphism, and conclude that 
either prolonged heating would yield a much greater decomposition of the magnetic 
constituent, or that a combination of pressure and heating is essential to produce 
the effects of the regional metamorphism. 

The Q value is somewhat lower than usual for this type of rock (eg. with the 
Tholeiites and Whin Sill of England which average 8 to 9) and this in part explains 
the low initial T value. The Z ratio averages 4.5 which indicates that the magnetic 
field in which the rock cooled in Pre-Torridonian times was at least 0.12 gauss, and 
that the rate of demagnetisation has been very slow, 80% in 800 million years. 

The S-2024 and S-100 samples have not changed susceptibility but increased 
their remanent magnetism by 2.0 and 1.7 respectively yielding Q values of 1.2 and 
1.0. The resultant Z ratios are 2.0 and 1.7 respectively, but since the remanent 
magnetisms are still very low these ratios may have a large percentage error. 

The S-243 samples have behaved rather differently, increasing their suscepti- 
bility by 50% and their remanent magnetism by a mean factor of 16. The Q values 
average 2.5 and, due to the very low T ratios, the Z value averages 11. It is notice- 
able that the Q value is similar to that for the unaltered dyke rock and it may be 
presumed that in some way the metamorphism in this case either completely des- 
troyed the remanent magnetism or left the rock in such a state that it rapidly — 
demagnetised. 


Summary 


By far the most striking result of these experiments is the discovery that Pre- 
Cambrian rocks can be found which possess remanent magnetism. ‘The field in which 
these rocks cooled was of the order of that existing today and not less than 0.12 
gauss. The rate of demagnetisation assuming that the field was 0.545 gauss, was 80% 
per 800 million years. The results for the metamorphosed dykes is in good agree- 
ment with the known geology and petrology, namely that the metamorphism was such 
as to destroy about 987 of the magnetite content. It would be interesting to know 
if this magnetite which remains was in any way different from the majority which 
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was destroyed, since in the study of the Tholeiite dykes it was found that at least 
one, and possibly two, types of magnetic constituent occur as about 2-3% of the total 
magnetite content, and that they are resistant to alteration and weathering. It is 
not necessarily suggested that the material is the same in all cases, but it is not 
easy to see why 1-2% of the magnetite should remain unless it was protected in 
some way, either by its inherent nature or some geological condition. 

Since all the metamorphosed dykes contain sphene, it may be concluded that, 
under the conditions of the heat-treatment in the laboratory, sphene does not decom- 
pose into magnetite, at least, in any case, not in these rock-types. 


Nomenclature 


The symbols used above are an extension of those derived by Koenigsberger 
(3). X is the ratio of susceptibility after heating (K,) through the Curie point to 
that found naturally (K,,). 

S is the ratio of remanent magnetism after (P) to that before heating (R). T 
is the ratio of permanent (R) to induced magnetism (K,H) of samples as found in 
the field. Similarly Q is the ratio of thermoremanent magnetism (P) to the suscep- 
tibility after heating (K,H). 

Z is the ratio of Q/T, which is also the S ratio corrected for changes in sus- 
ceptibility due to heat treatment. The units of K, R, P are 10-*cg.s. All other 
values are dimensionless ratios. H, the earth’s magnetic field in London, England 
was taken as 0.45 Oersted. 
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